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ABSTRACT 

We reinvestigate the relationship between the characteristic broad-line region size C^blr) an d the Balmer 
emission-line, X-ray, UV, and optical continuum luminosities. Our study makes use of the best available 
determinations of /?blr for a large number of active galactic nuclei (AGNs) from Peterson et al. Using their 

■ determinations of /?blr for a large sample of AGNs and two different regression methods, we investigate the 
robustness of our correlation results as a function of data sub-sample and regression technique. Though small 
systematic differences were found depending on the method of analysis, our results are generally consistent. 
Assuming a power-law relation /?blr x L a , we find the mean best-fitting a is about 0.67 ± 0.05 for the optical 

l-H ' continuum and the broad Hj3 luminosity, about 0.56 ± 0.05 for the UV continuum luminosity, and about 0.70 ± 

Ck 0.14 for the X -ray luminosity. We also find an intrinsic scatter of ~ 40% in these relations. The disagreement 

■ of our results with the theoretical expected slope of 0.5 indicates that the simple assumption of all AGNs having 
on average same ionization parameter, BLR density, column density, and ionizing spectral energy distribution, 
is not valid and there is likely some evolution of a few of these characteristics along the luminosity scale. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert — Quasars: general 
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1. INTRODUCTION 

There is mounting evidence that massive black holes reside 
in the centers of most or all massive galaxy bulges. Under- 
standing the demographics and properties of these black holes 
will hopefully clarify their roles in galaxy formation and evo- 
lution, the ionization of the intergalactic medium, and more. 
The masses of the central black holes in nearby non-active 
galaxies have been measured using stellar and gas kinematics 
of the central regions (Tremaine et al. 2002 and references 
therein). In active galactic nuclei (AGNs), the technique of 
reverberation mapping (Blandford & McKee 1982; Peterson 
1993; Netzer & Peterson 1997) has been used to measure the 
light-travel-time delay over which broad emission line flux 
responds to continuum luminosity variations, and to thus de- 
duce the characteristic size of the broad-line region (BLR) 
around the central, photoionizing source. By assuming the 
emission lines are broadened primarily by the virial gas mo- 
tions in the gravitational potential of the central object, the 
BLR size and the line width then give an estimate of the mass 
of the central object (Peterson & Wandel 1999). 

Wandel et al. (1999) compiled reverberation-based BLR 
sizes for 17 Seyfert-1 galaxies and derived a scaling rela- 
tion between BLR size and AGN luminosity. Kaspi et al. 
(2000) 6 measured reverberation-based BLR sizes for 17 high- 
luminosity AGNs from the Palomar-Green sample of quasars. 
Combining these measurements with those of the 17 Seyfert 
galaxies, Kaspi et al. (2000) found that the BLR size scales 
with AGN optical luminosity as Rblr L 010±om . Assum- 
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ing that this scaling relation is universal to AGNs of all lu- 
minosities and redshifts, several recent studies have used the 
relation to estimate the central masses in large AGN samples 
using a single-epoch measurement of the luminosity and the 
line width (e.g., Wang & Lu 2001; Woo & Urry 2002; Grupe 
& Mathur 2004). Since the reverberation studies above were 
based on the Balmer emission lines, generally broad Hj3, and 
on the optical AGN luminosity, single-epoch estimates for ob- 
jects at redshifts z > 0.6 have had to rely either on IR obser- 
vations (e.g., Shemmer et al. 2004) or on attempts to extend 
the optically based size-luminosity relation to UV luminosi- 
ties and UV broad emission lines (McLure & Jarvis 2002; 
Vestergaard 2002). While important progress has been made, 
there are still a number of potential problems that need to be 
addressed (e.g., Maoz 2002; Shemmer et al. 2004; Baskin & 
Laor 2005). 

In a recent study, Peterson et al. (2004) compiled all avail- 
able reverberation mapping data, obtained over the past 15 
years, and analyzed them in a uniform and self-consistent way 
to improve the estimates of the black-hole masses. They find 
that for a given luminosity, the black-hole mass can be pre- 
dicted to a precision (random component of the error) of typ- 
ically about 30%. However, there is also a systematic com- 
ponent of the error, due to the unknown geometry and kine- 
matics of the BLR, that amounts to a factor of about three, 
based on the scatter in the relationship between black hole 
mass and bulge velocity dispersion (Onken et al. 2004). In 
the present paper we use the BLR sizes compiled by Peterson 
et al. (2004) and study their relationship with the AGN lumi- 
nosity measured in X-ray, UV, and optical continuum bands 
and in the broad Hj3 emission line. In § 2 we explain how 
the data were derived, in § 3 we present the relation of the 
BLR size to the variously defined luminosities, and in § 4 we 
discuss our results. In § 5 we present a short summary. Un- 
less otherwise noted all wavelength mentioned in this work 
are common rest-frame wavelengths. 
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2. DATA 

Delay times for 35 AGNs with single or multiple data sets 
were calculated by Peterson et al (2004). For the current study 
we use the Balmer line time delays, as listed in their Table 6, 
Column 3 ("rest-frame time lags")- As noted by Kaspi et al. 
(2000), the derived BLR size from the first three Balmer lines 
(Ha, H/3, and Hy) are consistent with each other and aver- 
aging them reduces the scatter in the relation. We have cal- 
culated the error-weighted mean BLR size derived from these 
lines where available. We used only results which were desig- 
nated by Peterson et al. (2004) as reliable (see their section 4). 
While we do not necessarily expect each of the Balmer lines 
to have the same lag, we find empirically that the precision 
of the lag measurement for any one line is sufficiently low 
that averaging the Balmer line lags for any particular source 
reduces the scatter in the radius-luminosity relationship. We 
find a similar result when we average multiple H/3 lag and 
optical luminosity measurements for a single source. 

Note that in their final mass-luminosity relation Peterson et 
al. (2004) excluded several data sets for which the full width 
at half maximum (FWHM) determination was deemed to be 
unreliable. Since FWHM does not enter the present analysis 
we have included several of these data sets if their measured 
time lag was designated as reliable. We have also included an 
additional data set for NGC 405 1 presented by Shemmer et 
al. (2003). The mean BLR size for each data set is listed 
in Table [2 Column 3. The described analysis below was 
done twice - once using the Balmer-line averaged BLR size 
versus the luminosity, and a second time using only the H/3 
BLR sizes. The two analysis sets are discussed in § |3j In 
the case where only the H/3 BLR sizes were used, 3 objects 
(PG 0844+349, PG 121 1+143, and NGC4593) were excluded, 
since their H/3 -based BLR size measurement is considered 
unreliable. 

In the following subsections we explain how the luminosity 
at each wavelength band was calculated. We assume a stan- 
dard cosmology with Hq = 70 kms -1 Mpc~', Q.m = 0.3, and 
£l\ = 0.7. The luminosities are tabulated in Tables ^ an d |2] 
together with the measured fluxes (in the observer's frame) 
from which they were derived. 

2.1. General Considerations for the Multiwavelength 
Datasets 

We derived optical luminosities directly from the data sets 
from which the time lags were measured. Hence, the mea- 
surement of the time lag and the optical luminosity is si- 
multaneous. However, the UV and X-ray fluxes were mea- 
sured simultaneously with the BLR size for only a few AGNs. 
In these objects, we used these simultaneous UV and X- 
ray data in the analysis below (UV - NGC 3783, NGC 5548, 
NGC 4 1 5 1 , 3C390. 3, NGC7469, Fairall 9; X-ray - NGC4 151, 
NGC7469). Luminosities were not corrected for host galaxy 
contribution. To subtract this contribution, a knowledge of the 
spatial distribution of the host-galaxy luminosity is needed, 
combined with the slit width and seeing in which the obser- 
vation was done. This task is beyond the scope of the current 
paper. However, we note that this correction would be negli- 
gible for the quasars in the sample (AL^(5100A)> 10 44 ergs 
s _1 ) and it is likely small for Seyferts as well (e.g., Peterson 
et al. 1995; Kaspi et al. 1996). Since there might be a sys- 
tematic effect of greater host-galaxy luminosity contribution 
in lower luminosity objects, such a correction will tend to flat- 
ten the BLR-size - luminosity relation causing the power-law 



slope to be smaller. 

Recent studies (e.g., Gaskell et al. 2004; Hopkins et al. 
2004) indicate AGNs have reddening in addition to the Galac- 
tic reddening, probably due to gas in their host galaxy or the 
intergalactic medium. This is affecting the measured luminos- 
ity of the AGNs, though it is difficult to determine the intrinsic 
reddening of individual objects. Hence, we do not attempt to 
correct for such reddening in the study presented here. 

AGNs vary with differing amplitudes and timescales at dif- 
ferent wavelengths. Specifically, amplitudes become larger 
and timescales shorter as one progresses from the rest-frame 
optical, through the UV, to X-rays (e.g., Edelson et al. 1996). 
Therefore, another basic limitation of the comparison we 
make of the size-luminosity relation for luminosities in dif- 
ferent bands, beyond the fact that the actual sampling window 
in each band is very different (e.g., robust median luminosi- 
ties in the optical, based on tens of epochs, compared with 
few- or single-epoch-based luminosities in UV and X-rays), 
is the different statistical properties of the variability within 
each band. 

Table [fllists all multiple data sets. One line per object (in 
bold font) gives the mean of all the data sets for this object 
(cases in which there is only one data set appear with a normal 
font). This mean will be further discussed in §|3] 



2.2. Optical luminosities 

Optical luminosities were calculated from the mean 5100 A 
flux density of each data set. These luminosities are the same 
as the ones in Peterson et al. (2004; see their § 8), but with 
an improved correction for Galactic absorption. We used 
the E{B -V) values listed in the NASA/IPAC Extragalactic 
Database (NED) which are taken from Schlegel et al. (1998) 
and the extinction curve of Cardelli et al. (1989), adjusted to 
Ay /E(B — V) =3.1. For most objects the difference in lumi- 
nosities between the current work and Peterson et al. (2004) 
is less than 1%, while for a few it reaches up to 7%. The 
measured 5100 A flux densities and luminosities are listed in 
Table [2 Columns 6 and 7, respectively. The luminosity un- 
certainty is determined as the rms of the light curve of each 
data set. 

We have computed the H/3 luminosity from the average 
H/3 flux of each data set, correcting the fluxes for Galac- 
tic absorption as above. The H/3 flux and luminosity are 
listed in Table [2 Columns 4 and 5, respectively. We have 
also estimated the H/3 luminosity with the narrow compo- 
nent removed from the line profile (these narrow-line cor- 
rected H/3 fluxes and luminosities do not appear in TableQ. 
The contribution of the narrow component of H/3 was calcu- 
lated using the H/3 (narrow component )/[0 III] A5007 line ratio 
from Marziani et al. (2003), multiplied by the [O III]A5007 
line flux we measured in the data sets in hand. We used 
our own estimates for the narrow component of H/3 for: 
a) AGNs which are not in the Marziani et al. sample 
(NGC3227, NGC4151); b) AGNs which were observed at 
a low-state and for which we have a particularly good esti- 
mate of the narrow component of H/3 (Mrk79, PG 1229+204, 
NGC 5548, PG 1426+015). We also did not use the Marziani 
et al. estimates for five objects (PG 0804+761, PG 0844+349, 
PG 121 1+143, PG 1226+023, PG 1617+175) for which we 
think their [O III]A5007 line fluxes are overestimated due to 
blending with Fe IIA5018. The H/3 narrow component line 
flux we use in our analysis is listed in Table|2] Column 9. 
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Optical Observed Fluxes and Luminosities 
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0.386 ±0.060 




17 




7.41 ±0.95 


0.0523 ±0.0067 


8.12±0.91 


0.297 ±0.033 




17 


26.5 +4 - 3 


10.27 ±1.04 


0.0725 ± 0.0073 


13.47± 1.45 


0.492 ±0.053 




17 


24.81 3 ; 2 


9.34 ±0.61 


0.0659 ±0.0043 


11.83 ± 1.82 


0.432 ±0.066 




17 


6-511;? 


6.27 ±1.22 


0.0443 ±0.0086 


6.98 ±1.20 


0.255 ±0.044 
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TABLE 1 

— Continued 



Object 


Ref. a 


Time lag b 


F(Hj3) 

io- 13 


L(Hj3) 
10 43 


A(5100A) 
io- 15 


AL X (5KX)A) 
10 44 


(1) 


(2) 


lt-days 

(3) 


ergscm~ 2 s _1 
(4) 


ergs s -1 
(5) 


ergs cm -2 s~' A~' 
(6) 


ergs s _1 
(7) 




17 


14.3«j 


5.26± 1.12 


0.0372 ±0.0079 


7.03 ±0.86 


0.257 ±0.032 




27 


16.2+5.9 


7. 10 ±0.60 


0.0501 ±0.0042 
0.0532±0.0065 


5.78±0.62 


0.211 ±0.023 
0.334±0.041 


PG 1426+015 


2 


84 7+ 21 ' 4 


3.09 ±0.28 


0.6302 ±0.0581 


4.62±0.71 


5.21 ±0.80 


Mrk817 


1 


19.0+3* 


4.73 ±0.26 


0.1 100 ±0.0061 


6. 10 ±0.83 


0.75 ±0.10 




1 


153 -g 


4.00 ±0.40 


0.0930 ±0.0092 


5.00 ±0.49 


0.61 ±0.06 




1 


33.61" 
22.6+$ 


3.38±0.19 


0.0785 ±0.0045 


5.01 ±0.27 


0.612 ±0.033 








0.0938±0.0066 




0.656±0.065 


PG 1613+658 


2 


40 1+ 150 


2.02 ±0.1 3 


0.9496 ±0.0632 


3.49 ±0.43 


9.5±1.2 


PG 1617+175 


2 


86 - 8 -20l2 


1.40 ±0.17 


0.5 123 ±0.0636 


1.44 ±0.25 


3.00 ±0.52 


PG 1700+518 


2 


251 - 8 -38i8 


1.89 ±0.10 


5.5520 ±0.2882 


2.20 ±0.15 


42.34 ±2.89 


3C 390.3 


21 


23.61" 


2.12±0.17 


0.1950 ±0.0153 


1.73 ±0.28 


0.87 ±0.14 


Mrk509 


1 


79.6l«i 


12.30 ±0.06 


0.3985 ±0.0020 


10.92 ±1.99 


1.88 ±0.34 


PG 2 130+099 


2 


177 1+ 19 ' 9 
1 "- i -12.7 


5.12±0.46 


0.5558 ±0.0498 


4.84 ±0.45 


2.85 ±0.26 


NGC7469 


23 


4 5+06 
+' J -0.7 


7.80 ±0.35 


0.0577 ±0.0026 


13.57±0.61 


0.521 ±0.023 



NOTE. — Bold font lines give average values for objects with multiple data sets. The mean time lag is computed as a simple 
mean and its uncertainty is the rms of the individual measurements. The mean luminosity is computed as a simple mean and 
its uncertainty is the simple mean of the individual uncertainties (which are derived from the rms of each light curve). 



"Reference numbers are as in Table 1 of Peterson et al. (2004) and are given here to facilitate the identification of the data 
sets; reference 28 is Shemmer et al. (2003). 

b These are rest-frame time lags. For objects in which more than one Balmer-line time lag was measured we list the averaged 
rest-frame time lag. 
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TABLE 2 
Object Characteristics 



UDJcCl 


D ar\ c lii ft 

Keasnin 


77 ( J} \/\ 
C, [ D — V ) 




3 7" 1 1 A SO A \ 


r (^Z— 1U Kc v ^ 


t so in 
L\L— iu Ke V ^ 


Kcr. 


17/11 ft 








10 


1 n44 
10 


10 


1 a43 

10 




10 








—9 — 1 S — 1 

ergs cm s A 


_ i 

ergs s 


—9 — 1 

ergs cm s 


_ 1 

ergs s 




ergs cm * 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


MrV IIS 


09S78 
U.UZJ / o 


035 


7 2 ± 1 


9 1 9 -1- 10 

Z. 1Z x U.jU 


Q 41 4- 10 

y.H-J ZTZ U. 1U 


1 4.in 4- 01 S 
i ,4ju zrz u.u i j 




410 


pn 0096+1 9Q 

L V 1 UUZ.UT 1 Z.? 


14900 

VJ. I'rZ.UU 


071 


1 5 ±0 5 


97 i£ _l 7 4.5 


u. jj ziz v J . i J 


1.A 1 6 4- 70 
ji. i u zrz vj. / vj 


3 


063 


pn 00S9+9S1 


1 SSOO 


047 


1.1 ±0.2 


1 6 88 -4- 1 07 
iu.oo zrz j .vj / 


1 OQ 4- 17 

i v ) . v / y zrz u.i/ 


66 ± 1 1 




091 


1 all till 7 


04702 


027 


3.5 ±1.0 


a qo _i_n 07 
j.jo zrz vj. y i 


1 Q Oft 4- 1 00 

17.UUZTZ 1 .UU 


q 09 _4 a 59 
7.7Z zrz v j . j z. 


9 


283 


Mrle SQO 

IVllJv J;?U 


09618 

V7.UZ.UJO 


037 


2 4± 1 4 


75 _|_n A A 
u . i j zrz vj . tt 


97 00 4- 9 70 
_ / . uu zrz z. , / u 


4 29 ± 43 


7 


121 


3C 120 


03301 


297 


1.2± 1.2 


4 45 ±4 23 


4.4 ?8 4-9 99 
ii . z,o zrz z..z.z. 


11 1 1 4-n S6 

1 1 . i j zrz vj . j u 


2 


151 


Aim 190 

rtJVll 1 ZAJ 


01910 


128 


3 5 ±0 5 


a _i_ n as 
j.jt zrz vj.*to 


97 14. -U 9 71 
z/.ji zrz z.. / j 


£ 57 4. a 66 
u.j / zrz u.uu 


5 


364 


Mrk79 


02219 


0.071 


2.0± 1.0 


57 + 29 


9< q<; + 9 f,0 

.y\j i z.uvi 


2 Q(i + 29 


7 


0.221 


pn 0804-u76i 


10000 
VJ . 1 uuuu 


035 


7 ± 1 


jo.jozr j. z.z 


8 14 4- 07 
o . jh- zrz u.u / 


91 114-018 
Z.1.U zrzu.io 




VJ. 


pn 0844+14Q 

1 V 1 UOTTTJ77 


06400 


037 


1.7 ±0.5 


5 114-1 00 

J .~t 1 _i_ 1 .UU 


9 37 4. n 09 
z.. j / zrz u.uz. 


9 14. 4- o 090 
z,. j*t zrz vj.uz.vj 






Mrk- 1 1 


01S9Q 

U.U J JjL7 


013 


1 0±0 5 


48 4-0 94 

U.tO ZTZ U.Z.'t 


96 64 4- 1Q 

Z. U . Irr ZTZ u.17 


7 fi7^ -J_ n OSS 
/ .vj / u zrz vj.uj j 




148 


PG0953+414 


0.23410 


0.013 


2.0± 1.0 


64.3 ± 32.1 


2 5» + o 04 


42 22 + 65 


2 


0.534 


NfiC 3227 


00386 


0.023 


n 04 ±0 01 

\j . ut i V7.ui 


0003 + 000 1 


25 52 + 1 30 

AJ.JZ 1 1 . J U 


0839 + 0043 

U.UOJ7 — J— U.UU^J 


2 


0.558 


wr.r i< i ft 


00884 


042 


3 0± 1 5 


n 10ft 4- 0S1 

U. 1UO ZTZ U.UjJ 


7Q A8 -U 1 70 
/ j .vjo zrz j . / u 


1 980 4- 064 
i .z.oU zrz U-UO^t 


2 





Krnn 1781 


00Q71 
u.uu:j / J 


119 


4g + 10 


n 171 4- 078 
u.j / j zrz u.u / o 


J 1 .*-ry ZTZ z.JJ 


i 08S 4- 0S4 
i .uoj zrz u.uj^t 


2 


656 


Nnn 40s i 


00714 

U.UUZ. J*"T 


013 


1 2±0 2 


n 0091 4-0 0004 

U.UUZ.J ZTZ U.UUU'T 


21.21 ± 1.10 


096 4-0 0011 
u . uz.u zrz u . uu i j 


2 


o 


Nnn 41 si 

1> V I v . '-t 1 J 1 


o 0011 9 

VJ.UU J JZ. 


028 


35. ± 15. 


n 1 <A _|_ n n^A 

U. 1 JT ZTZ U.UUU 


201. ± 10. 


4Q0 4- 09 S 
vj .^y\j zrz u.uz.j 


2 


9 912 


pn 1?1 1 -1-1 41 

.TVj IZI I+Itj 


080Q0 
U.UoUVU 


035 


3 ± 1 


Q O _L T Q 
V.O Zl J.J 


9 89 4- 04 
z.oz zrz u.u*+ 


4 SS7 4- 06S 
t.jj / zrz u.uoj 




Q 


Pn 1 996+091 


15834 


021 


15 + 40 


900 6 4- SI S 

ZW.U ZTZ JJ.J 


qn 11 -4 90 8 
?vj. jj zrz z.u.o 


619 ± 142 


3 


o 


Pn 1 99Q+904 


06101 
vj.uu ju i 


027 


2.0± 1.0 


3.6 ± 1 8 


9 54 -4 16 
i.ji zrz v / . i u 


2.43 ±0.15 




0.152 


Nnr 4sqi 


00Q00 


025 


2.0 ± 1.0 


o&a a. n 019 

U.Ulrr ZTZ U.UJZ 


11 RQ 4- 1 60 
j i .oy zrz i .uu 


n 575 4. n 09Q 
u. J / J zrz u.uz.? 


2 


217 


pn i 107+08S 

tKJ IjU/tUoj 


1 SSOO 
U. 1 JJUU 


014 


8 + 02 


11 2 ± 2 8 


1 6S 4- IS 
j .vjj zrz u.j J 


23 9 ± 2 3 




155 


TP 419QA 


01 60S 
U.U 1 OUJ 


OSQ 
U.UJ y 






77 Af. _j_ i on 
/ / .'+o zrz j.7U 


4. 4.Q _4 n 91 
t +. L ty! zrz U.z. J 


2 


253 


MrV 970 


0104 S 

VJ.VJ-.lVyT J 


016 


6 + 05 


99 _i_n 1 Q 

U.Z.Z. ZTZ U. 1 O 


95 09 _4 9 5Q 

ZJ.7Z ZTZ Z.J7 


5 59 _4 o 55 

j . j z. zrz v j . j j 


5 


325 


PG 1411+442 


0.08960 


0.008 


1.4 ±0.4 


4.70 ±1.34 


0.78 ±0.09 


1.57±0.18 


1 


0.140 


NGC 5548 


0.01717 


0.020 


3.0±1.0 


0.34 ±0.11 


46.05 ±2.30 


3.06±0.15 


2 


0.614 


PG 1426+015 


0.08647 


0.032 


7.5±2.0 


27.82 ±7.42 


5.25 ±0.25 


9.77 ±0.47 


3 


0.072 


Mrk817 


0.03145 


0.007 


1.5 ±1.0 


0.54 ±0.36 








0.078 


PG 1613+658 


0.12900 


0.027 


3.5±1.0 


30.52 ±8.72 


6. 14 ±0.09 


26.90 ±0.39 


3 


0.041 


PG 1617+175 


0.11244 


0.042 


0.75 ±0.35 


5.35 ±2.50 








0. 


PG 1700+518 


0.29200 


0.035 


0.6 ±0.4 


39.4 ±26.3 








0. 


3C 390.3 


0.05610 


0.071 


0.6±0.2 


1.17±0.39 


22.80 ±5.00 


17.11 ±3.75 


4 


0.172 


Mrk509 


0.03440 


0.057 


5.5±1.0 


3.48 ±0.63 


43.10±2.10 


11.79±0.57 


2 


0.777 


PG 2 130+099 


0.06298 


0.044 


2.2±0.2 


4.50 ±0.41 


4.80±0.12 


4.59±0.11 


3 


0.385 


NGC 7469 


0.01632 


0.069 


3.0±0.5 


0.451 ±0.075 


27.64 ±3.40 


1.66 ±0.20 


8 


1.357 



REFERENCES. — FOR COLUMN 8:1- George et al. (2000); 2 - George et al. (1998); 3 - Lawson & Turner (1997); 4 - Leighly et al. (1997); 5 - 
The Tartarus Database lhttp://tartarus. gsfc.nasa.gov/ 1; 6 - Williams et al. (1992); 7 - Turner & Pounds (1989); 8 - Nandra et al. (1998). 

"Observed X-ray flux in the observed energy band 2-10 keV. 
b X-ray luminosity in the rest-frame energy band 2-10 keV. 



2.3. UV luminosity 
For the UV luminosity, we have used both the 1350 A 
(used by Vestergaard 2002) and 1450 A (used by Shemmer 
et al. 2004) wavelength ranges. UV fluxes at these rest frame 
wavelengths were measured, for most objects, using FOS and 
STIS spectra taken from the Hubble Space Telescope (HST) 
archive. In cases where multiple observations exist, we have 
used the averaged flux density, except for the few cases where 
the UV data were obtained simultaneously with the optical 
time lag measurements - see §|2] In several cases (Mrk 110, 
PG 1229+204, PG 1426+015, Mrk817, PG 1613+658) there 
were no observations at 1350 A or 1450 A but spectra at about 
1200 A could be extrapolated (using a straight line) in order to 
estimate the 1350 A and 1450 A flux densities. For several ob- 
jects without HST UV spectra we used archival International 
Ultraviolet Explorer (IUE) spectra to measure the 1350 A and 
1450 A flux densities. These objects are 3C 120, Mrk590, 
Mrk 79, and PG 1617+175. For one object, IC4329A, no use- 
ful UV data were found. The flux densities were corrected for 
Galactic absorption, as described above for the 5100 A flux 
density. The measured 1450 A flux densities and calculated 



luminosities are listed in Table|2|in Columns 4 and 5, respec- 
tively. 

2.4. X-ray luminosity 

Most of the objects with reverberation mapping data also 
have X-ray data in the literature and in archives. We searched 
the literature for the 2-10 keV fluxes of each object (see ref- 
erences in Table|2j column 8). For 3 objects (PG 1700+519, 
Mrk817, and PG 1617+175), no useful X-ray data were 
found. The fluxes were corrected for Galactic absorption us- 
ing the PIMMS tool (Mukai 1993) and the published power- 
law slopes and Galactic column densities. The observed 2- 
10 keV band was K-corrected to yield the rest-frame 2-10 
keV luminosity. We note, however, that because of the small 
Galactic absorption of most objects and their small redshifts, 
these corrections are small (of order of a few percent and 
hence much smaller than the known X-ray variability). Also, 
any intrinsic absorption has a negligible effect in this hard 
X-ray energy range. The measured fluxes in the observed 2- 
10 keV band and the calculated luminosities in the rest-frame 
2-10 keV band are listed in Table|2]Columns 6 and 7, respec- 
tively. 
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3. BLR SIZE - LUMINOSITY RELATIONS 

The relations between the BLR size and the various lumi- 
nosities described above are plotted in Figure ^ The power- 
law relation between the BLR size and XL^ (5 1 00 A) is clearly 
confirmed by our analysis. Such a relation is also apparent 
when UV, X-ray, and H/3 luminosities are used. Two obvi- 
ous outliers, which do not fit the general trends in any of the 
wavelength bands, are NGC3227 and NGC4051. These are 
the two-least luminous objects in all wavelength bands. This 
behavior could be due to major reddening or obscuration in 
these two objects, which decreases their measured luminosity 
(e.g., NGC 3227 - Crenshaw et al. 2001; NGC 4051 - Kaspi et 
al. 2004). For NGC 3227, the extinction correction suggested 
by Crenshaw et al. (2001) increases the UV luminosity (at 
1450 A) by a factor of ~ 66 and the optical luminosity by a 
factor of ~ 2 (see Figure^. Although this would improve the 
agreement of this data point with the general trend, it would 
still remain an outlier. An alternative explanation can be that 
the BLR-size - luminosity relation found at high luminosi- 
ties breaks down at low luminosities. We will exclude these 
two objects from our fits to the BLR-size - luminosity rela- 
tion, noting the derived relation applies to AGNs with optical 
luminosity in the range 10 43 < AL A (5100A) < 10 46 ergs s _1 . 

We have used two methods to calculate the relations be- 
tween the BLR size and the various luminosities: 

1. The linear regression method of Press et al. (1992), in 

which a straight-line is fitted to the data with errors in 
both coordinates (known as FITEXY). This method is 
based on an iterative process to minimize % 2 . We fol- 
low Tremaine et al. (2002), who account for the intrin- 
sic scatter in the relation by increasing the uncertain- 
ties 7 until obtaining % 2 per degree of freedom equal to 
1 (see below in Table[3]column 5). 

2. The bivariate correlated errors and intrinsic scatter 
(BCES) regression method of Akritas & Bershady 
(1996). This method takes into account the uncer- 
tainty in both coordinates as well as the intrinsic scatter 
around a straight line. In our analysis, we use the BCES 
bisector result which is the mean of the two fits: one 
when fitting Y as dependent on X, and the other when 
fitting X as dependent on Y. 

Tremaine et al. (2002) compare the two methods using 
Monte-Carlo simulations and concluded that the FITEXY is 
superior to the BCES method. However, in their implementa- 
tion for the BCES method Tremaine et al. (2002) did not use 
the BCES bisector value. They only used the fit of Y as depen- 
dent on X. This might have distorted their conclusion when 
comparing the two methods. Thus in the analysis below we 
use both methods. One advantage of the BCES method over 
the simple FITEXY method is its ability to account for real 
intrinsic scatter, if it exists in the data. However, the available 
code of the BCES method does not give any quantification 
of this scatter. The implementation of the FITEXY method 
that we use (of Tremaine et al. 2002) also considers intrinsic 
scatter and allows us to quantify it. Such intrinsic scatter was 
assumed in previous investigations of the BLR-size - lumi- 
nosity relation (e.g., Vestergaard 2002), but not quantified. 

The two fitting methods do not account for asymmetric un- 
certainties (e.g., as we have for the uncertainties of the BLR 

7 In the present work we increase the uncertainty of the time lags by a 
certain percent of the measured value. 



sizes). In order to account for the asymmetric uncertainties, 
we used in our fits the uncertainty value in the direction of the 
fitted line. This typically required a few iterations to converge 
completely. 

The correlation analysis below is done in two ways: (1) by 
using all data sets that are available, where in case of mul- 
tiple data sets for the same object we treat each data set as 
independent; (2) by computing the mean BLR size and mean 
luminosity of all the multiple data sets of the same object, and 
using only that mean for each object. 

In Table [5] we summarize the results from the different fits 
we carried out for the relation of the various luminosities with 
the BLR size. For each of the luminosities measured, two 
lines are given: the first line lists the results when using all 
data sets and the second line lists the results when using the 
averaged data set per object. The number of data sets used 
in the fit is listed in Column 2. The relations are given in the 
form: 

Rblr b m 

101tdays~ AL (i) 

where L is the luminosity normalized as given in column 1 of 
Table|3] The parameters A and B are given in columns 3 and 4 
for the FITEXY method and in columns 6 and 7 for the BCES 
method. In column 5 we list the intrinsic scatter, in percent 
of the value, found when using the FITEXY method. Also, in 
columns 8 and 9 we list the Pearson and Spearman correlation 
coefficients, respectively. 

As explained in §|2] the analysis was done twice - once us- 
ing the mean Balmer-line BLR size versus the luminosity (re- 
sults are presented in the upper part of Table|3jl, and a second 
time using only the H/3 BLR sizes (lower part of Table[3}. The 
results from the two analysis sets are consistent and no signif- 
icant differences are detected. The correlations coefficients 
(columns 8 and 9 of Table [5} are somewhat higher when us- 
ing the mean Balmer-line BLR size and this may indicate that 
using these lags reduces the scatter in the relation with the 
luminosity. Thus, in the analysis below we refer only to the 
upper part of Table[3]which uses the mean Balmer-lines time 
lags. 

The results of Table[3]indicate that there is a systematic dif- 
ference between the two fitting methods. The BCES method 
slope parameter (B) is generally larger than the one found 
with the FITEXY fit, and also the constant parameter (A) is 
generally larger. This illustrates the important issue of how 
the regression method may influence the correlation results; 
using different regression methods yields different regression 
slopes. We could not find what is the main cause for the differ- 
ences and we attribute them to the different algorithms, strate- 
gies, and statistical concepts each of the two methods uses. 
The real uncertainty of the slope measurement should be the 
range between the different methods and not necessarily the 
formal error given by each method. Furthermore, the regres- 
sion methods are sensitive also to outlier points. We find that 
the BCES method is much more sensitive to the inclusion or 
exclusion of the two lower-luminosity points (NGC 3227 and 
NGC 405 1 ) compared to FITEXY. 

In its present implementation the FITEXY method estab- 
lishes that intrinsic scatter does exist in the BLR size - lumi- 
nosity relation. Only when we artificially increase the formal 
measured uncertainty in the BLR size do we find a % 2 of unity. 
We find that the intrinsic scatter in the relation is ~ 40% (Ta- 
ble [3] Column 5). We note that this intrinsic scatter can be 
caused by a real scatter in the physical BLR size - luminosity 
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FlG. 1 . — Balmer-line averaged BLR size plotted versus various luminosities given in ergs s -1 . The left side panels show all data sets and the right side panels show one averaged 
data set per object. The data are listed in TablesfTlandpl For each panel two fits (excluding the two lowest luminosity objects) are shown: solid line - FITEXY and dashed line - BCES. 
Fits parameters are given in the upper part of Table|3| Right pointing arrows show the shift in location of the data points of NGC 3227 when the Crenshaw et al. (2001) extinction 
correction is applied. 
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relation but also by some systematic effects such as intrinsic 
reddening, contributions by the host galaxy, or effects of vari- 
ability due to non-contemporaneous observations. There is 
currently no way to distinguish between the different contrib- 
utors to this scatter. 

In Figure 13 we plot the mean Balmer-line BLR size ver- 
sus the AL^(5100 A) luminosity, with one averaged point per 
object. We also plot four fits: using all 35 points, excluding 
the low luminosity AGNs, and with the two fitting methods. 
Within the luminosity range of the measurements (10 43 — 10 46 
ergs s _1 ) all fits are consistent with each other and all are 
well within the scatter of the points in the plot. Also, when 
using only the data for the more luminous (>5.2x 10 43 ) ob- 
jects (26 objects) from Figure|2] we find a slope of Bfitexy = 
0.63 ±0.07 and a constant of Afitexy = 2.4 ±0.3, which are 
consistent with the results when using the whole sample. 

Comparing the two methods of analysis, using all data sets 
versus using only one averaged data set per object, we find 
that the slope parameter (B) is somewhat different while the 
normalization parameter (A) is similar (both for the BCES 
and FITEXY methods). In general, when using one averaged 
data point per object, the slopes are larger by ~ 0.05 (which is 
within the formal uncertainty) than when using all data sets. 
An argument for using one averaged data point per object is 
that we do not know if the BLR-size - luminosity relation in 
a given object is the same as between different objects (the 
physical processes driving the relation in a single AGN might 
be different from the ones driving it in a range of objects). 
Thus, in order to have a relation that spans a large range of lu- 
minosity and different objects and to study the physics across 
the luminosity function, it is better to use the relation with a 
mean BLR-size and a mean luminosity for each object. On the 
other hand, if the relation in one object is the same as between 
objects there is justification to use all data sets. This has been 
tested observationally for only one object, NGC 5548, by Pe- 
terson et al. (2002) who found the BLR-size - luminosity 
relation has a slope of 0.95 when using optical luminosity and 
0.53 for the UV luminosity. 8 While the slope with the UV 
luminosity is similar to the one we find in Table|3] the optical 
luminosity slope is not. Hence, it is still not clear if the re- 
lation in a single object is the same as in a sample of objects 
and we present both results in the current study. 

Removal of the Hj3 narrow component from the total H/3 
luminosity affects the fit results only at the 1% level with 
no obvious general trend. Also, when correlating the BLR 
size with the UV luminosity, the difference between using the 
1350 A or the 1450 A bands is negligible (less than 2% and 
well within the uncertainties), as expected. The flux densities 
in these bands are very similar, and in some cases identical. 

Wu et al. (2004) correlated the BLR size with the H/3 lu- 
minosity for the objects of Kaspi et al. (2000) using the data 
given there. The H/3 luminosity can be considered a good in- 
dicator for the AGN luminosity as the line is driven by the 
ionizing continuum. It was noted by Kaspi et al. (2000) that 
the use of the H/3 luminosity instead of the optical luminosity 
do not change their correlation results (see their §4.2). Wu et 
al. (2004) use only one data set per object and assume a dif- 

8 The 0.53 slope was found by scaling the relation for the optical lumi- 
nosity (0.95) by the slope of the relation between the UV and the optical 
luminosities, in which Peterson et al. (2002) find to be 0.56. We note, how- 
ever, that Gilbert & Peterson (2003) favored a slope of 0.67 for the relation 
between the UV and the optical luminosities, in which case the BLR-size - 
UV luminosity slope of Peterson et al. (2002) is then 0.64. 



ferent cosmology (Hq = 75 km s _1 Mpc -1 and qo ~ 0.5) 9 and 
obtain a slope of 0.68 ±0.11. This result is consistent with 
the one we obtain in our analysis (Table[3}. 

4. WHAT DETERMINES THE BLR SIZE? 

A theoretical slope for the relation between the BLR size 
and the luminosity, predicted by various early papers, is 0.5. 
This is based on the assumption that, on average, all AGNs 
have the same ionization parameter, BLR density, column 
density, and ionizing spectral energy distribution (SED). A 
slope of 0.5 is also expected if the BLR size is determined 
by the dust sublimation radius (Netzer & Laor 1993). How- 
ever, several studies have indicated that the SED is luminosity 
dependent, and thus at least part of the assumptions above 
are not justified (e.g., Mushotzky & Wandel 1989; Zheng & 
Malkan 1993; Puchnarewicz et al. 1996). Here we use the ob- 
served relations of the BLR size with the different luminosity 
bands to study which band influences the BLR size the most. 

4. 1 . UV vs. Optical Continuum Luminosity 

Use of the H/3 and AL^(5100A) luminosities in the BLR- 
size - luminosity relation yields similar slopes of about 0.69. 
The slope with the UV luminosity [both AL^(1450A) and 
AL^(1350A)] is lower by about 0.1 than the slopes using the 
optical bands (this is manifested both in the FITEXY method 
and the BCES method). The UV-based slope is the closest 
to the theoretical expectation of 0.5 based on the assumptions 
mentioned above. 

The UV band (i.e., 1200-3000 A) is close in wavelength 
to the ionizing continuum which has a large effect on the re- 
sponding BLR gas. However, the UV luminosities are based 
on measurements that are generally not simultaneous with the 
reverberation measurements, and sometimes greatly separated 
in time. The luminosity in the optical band (i.e., 4000-8000 
A) might suffer from systematics which skew the BLR-size 
- luminosity relation. These systematics may be the stellar 
contribution to the optical luminosity of the Seyfert galaxies 
(which was not taken into account in this work), or a change 
of the AGN SED with luminosity, which translates to differ- 
ent slopes in the BLR-size - luminosity relation. The change 
in SED is suggested by the fact that brighter AGNs have flat- 
ter optical-UV continua, ranging from an average slope (/3 in 
L v «: v$) of about —1 for Seyferts, to about —0.3 for quasars 
(e.g., Neugebaueret al. 1989; Schmidt & Green 1983; Francis 
et al. 1991). This means that quasars have a larger UV/optical 
flux ratio than Seyferts (a quasar which is 10 times more lu- 
minous than a Seyfert in the optical band will be more lumi- 
nous than the same Seyfert by a larger factor in the UV band 
). Thus, when using UV, rather than optical luminosities, the 
ratio of luminosities between quasars and Seyferts grows and 
the slope of the BLR size versus luminosity relation becomes 
shallower as indicated by Table [3] 

Nonetheless, as already discussed and as we see below, po- 
tentially better indicators of the ionizing continuum, such as 
the H/3 luminosity, seem to contradict the finding that 0.5 
is the "correct" slope. For comparison we plot in Figure |2j 
showing the relation between mean BLR size and optical lu- 
minosity, a FITEXY fit to 33 data points with a slope fixed to 

9 This different cosmology would cause a slightly flatter slope for the BLR- 
size - H/3 -luminosity relation compared with the cosmology used in the cur- 
rent work by about 0.05, which is well within the uncertainty range of the Wu 
et al. (2004) result. 
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TABLE 3 
BLR size — Luminosity relation" 



Luminosity N Afitexy Bfitexy Scatter Abces Bbces Pearson Spearman 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



Using mean Balmer-lines time lag 



L(Hj3)/10 43 


59 


7 on+0.87 

'•"-0.78 


0.619±0.045 


38 


o 42+0-93 
y -^ z -0.85 


0.670 ±0.050 


0.885 


0.860 




33 


7 cc+0.90 
'• JJ -0.81 


0.687 ±0.058 


40 


8 91+ ' 92 
°' yl -0.83 


0.690 ±0.068 


0.916 


0.933 


L(Hj3)/10 43 no narrow 


59 


8 14+ ' 91 
°- 1# -0.82 


0.618 ±0.044 


38 


9 68 +0 59 
'■ D °-0.90 


0.659 ±0.052 


0.879 


0.872 




33 


7 07+0.89 
'■ S '-0.80 


0.683 ±0.055 


40 


Q 2S+0-99 
'• ZJ -0.89 


0.686 ±0.064 


0.916 


0.933 


Az^iooyio 44 


59 


24 5+o.l8 
z - 4;) -0.17 


0.608 ±0.045 


40 


2 65+0+9 


0.694 ±0.056 


0.869 


0.845 




33 


2 n7 +0.26 
z - u '-0.23 


0.665 ±0.065 


46 


2 on+0.32 
z - Jy -0.28 


0.723 ±0.075 


0.884 


0.910 


XL^(U50)I\0 U 


58 


2 12+0.17 
z - lz -0.15 


0.496 ±0.042 


41 


2 27+°' 17 
z ' z '-0.16 


0.595 ±0.046 


0.865 


0.807 




32 




0.545 ±0.063 


45 


2 ,+0.24 
z - ul -0.21 


0.604 ±0.053 


0.888 


0.910 


AL^nsoyio 44 


58 




0.499 ±0.040 


41 


7 7C+0.17 
Z - ZJ -0.16 


0.595 ±0.046 


0.862 


0.803 




32 


i 7C+O.26 
1 -' J -0.22 


0.540 ±0.056 


46 


2 07+0- 23 
z - u '-0.21 


0.584 ±0.046 


0.887 


0.897 


L(2-10 keV)/10 43 


54 




0.614 ±0.064 


52 


65+ 013 


0.830±0.102 


0.705 


0.669 




30 


0.86** 


0.544 ±0.091 


64 


7,+0.24 


0.813±0.119 


0.688 


0.701 



Using only H/3 time lag 



L(HjS)/10 43 


55 


7 55+ - 91 
'■ JJ -0.82 


0.609 ±0.048 


38 


47+1-02 
y '^ z -0.92 


0.669 ±0.056 


0.874 


0.845 




30 


7 1 s+0.94 
'• 1S -0.83 


0.694 ±0.064 


40 


8 79+ 100 
°' /y -0.90 


0.691 ±0.080 


0.906 


0.923 


L(Hj3)/10 43 no narrow 


55 


7 o S +0.94 
'■°-'-0.84 


0.604 ±0.047 


38 


n 73+I.IO 
'•'•'-0.99 


0.657 ±0.059 


0.867 


0.854 




30 


7 4Q+0.99 
'• Hy -0.88 


0.687 ±0.063 


40 


9 18+ 109 

y -'°-0.98 


0.685 ±0.081 


0.906 


0.925 


AL^siooyio 44 


55 


243+0.I8 
z - 4 - 5 -0.17 


0.601 ±0.047 


39 


2 69 +0 - 20 
Z - DV -0.18 


0.670 ±0.060 


0.870 


0.843 




30 


z - uu -0.24 


0.665 ±0.069 


44 


2 47+0-36 
z '^'-0.32 


0.674 ±0.089 


0.882 


0.911 


Az^cMsoyio 44 


54 


2 ,,+0.17 
z - 11 -0.16 


0.498 ±0.044 


39 


7 7S+0.17 
Z - ZS -0.16 


0.566 ±0.047 


0.858 


0.780 




29 


1 77+0.29 
1 -' z -0.24 


0.548 ±0.061 


44 


7 13+0.29 
z - 1 '-'-0.26 


0.552 ±0.062 


0.878 


0.894 


Az^nsoyio 44 


54 


2 06+o.is 

Z - UD -0.16 


0.495 ±0.041 


40 


2 26 +017 
z ' zo -0.16 


0.557 ±0.045 


0.854 


0.776 




29 


1 77+0.27 
1 -' z -0.23 


0.536 ±0.058 


44 


2 M+0.27 
z -' -0.24 


0.548 ±0.050 


0.876 


0.878 


L(2-10 keV)/10 43 


50 




0.612 ±0.063 


48 


65+ 015 


0.806 ±0.108 


0.700 


0.690 




27 


o-«ds 


0.532 ±0.090 


62 


u - ' !-0.22 


0.773 ±0.139 


0.668 


0.695 



NOTE. — For each luminosity in column 1 two rows are given: first row is the fit results where multiple data sets were used for each object, 
and the second row is the fit results where all data sets per object were averaged. See text for details. 

"Results of the fit which correspond to equation 1 . 
b Given in percent of the value. 



0.5. a line with a slope. The best fit A parameter (see equa- 
tion^ is 2.23. 

We note in passing, that for individual AGNs the ratio of 
the UV luminosity to the optical luminosity is observed to 
increase as the object becomes brighter (e.g., Peterson et al. 
2002; Gilbert & Peterson 2003). This can affect the BLR-size 
- luminosity relation for those AGN with multiple data sets. 
A possible prediction would be that the flattest slopes might 
result from BLR-size - luminosity relations that were derived 
from the case when including all individual data sets, rather 
than using the average for each AGN. 

4.2. Z/j3 luminosity 

Simple arguments based on recombination theory suggest 
that, in low-density low-optical depth gas the Hj3 emission- 
line luminosity is a good indicator of the ionizing continuum 
luminosity as the line luminosity is driven by this continuum. 
The physical conditions in the BLR are, however, known to 
be more complex (Netzer 1990 and references therein) yet we 
expect a general linear dependence between the ionizing flux 
and the Hj3 luminosity. The slope of the BLR-size - H/3 lu- 
minosity relation found here is almost identical to the slope 
found using AL A (5100A). This is in accord with studies that 
show the Balmer-line luminosity scales as the optical lumi- 



nosity between different AGNs (e.g., Yee 1980; Shuder 1981). 
However the BLR-size - H/3 luminosity relation does not 
agree with the slope obtained using the UV continuum. This 
seems to suggest that the A 1450 continuum is not the best in- 
dicator of the ionizing continuum or, perhaps, that our under- 
standing of the BLR physics is incomplete, i.e., that the H/3 
line luminosity may not depend linearly on the ionizing con- 
tinuum due to optical depth and/or collisional de-excitation 
effects in hydrogen lines. It has also been argued that a major 
line driver, and a possible BLR-size regulator, is the very high 
energy (X-ray) continuum. We do not favor this explanation 
because the fractional energy in this continuum is small com- 
pared with other wavelength bands. This is all supplemented 
by the fact that the observed line response reflects changes 
in ionization over a large volume and is geometry dependent 
in a way that obscures the local, immediate line response to 
continuum variations. For example, the covering factor and 
mean gas density of the BLR gas might decrease with distance 
from the central ionizing source (e.g., Kaspi & Netzer 1999), 
and the non-linear response in the inner BLR thus might con- 
tribute more than the response from the outer parts. Evidently, 
the BLR-size - luminosity relationship is rather complicated 
and depends on several factors. 
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^(5100 A) 

FIG. 2. — Balmer-line averaged BLR size plotted versus the AL^(51()0 A) luminosity (in units of ergs s~') on a larger scale. The BLR size of each data set is 
determined from the averaged Balmer-line time lags. Objects with multiple data sets have been averaged to one point per object. Five line fits are shown: dotted 
line - FITEXY using all 35 points; dot-dashed line - BCES using all 35 points; solid line - FITEXY excluding the 2 lower luminosity points; dashed line - BCES 
excluding the 2 lower luminosity points; dot-dot-dashed line - FITEXY excluding the 2 lower luminosity points and fixing the slope to 0.5. 

for the flux ratios between the X-ray and optical bands) anti- 
correlates with the luminosity (e.g., Vignali et al. 2003). This 
implies that when the optical luminosity increases the X-ray 
luminosity increases by a smaller factor, and this can explain 
the larger slope we see in the relation of the BLR size with the 
X-ray luminosity. 

5. SUMMARY 

Peterson et al. (2004) recently compiled and analyzed in a 
uniform and self-consistent way the BLR size from all avail- 
able AGN reverberation mapping data obtained over the past 
15 years. In this paper, we correlate the BLR size with the 
luminosity L, in the X-ray, UV, and optical bands. We in- 
vestigate the correlation with two different regression meth- 
ods as well as using multiple data sets for individual objects 
versus only one averaged data set for each object. Though 
small systematic differences exist between the different meth- 
ods of analysis, the results are generally consistent. Assuming 
a power-law relation Rblr ^ L a , we find the mean best-fitting 
a is about 0.67 ±0.05 for the optical continuum and the broad 
H/3 luminosity, about 0.56 ±0.05 for the UV continuum lumi- 
nosity, and about 0.70 ± 0. 14 for the X-ray luminosity. These 
slopes are a simple average of the slopes found from the dif- 
ferent methods using mean Balmer-lines time lags and the un- 



4.3. X-ray Continuum 

The BLR size versus X-ray luminosity relation shows a 
large scatter (50-60%), though a general trend can be found. 
The two regression methods give somewhat different results 
for the slope. We are not able to explain the reason for this 
difference though it is likely that it is caused due to the large 
scatter and the way each method utilize the data. The slope 
of the best fit relation is similar in the FITEXY method to 
the slope from using the optical luminosity but in the BCES 
method it is steeper by about 0.1-0.2 (the best fit slope is of 
order ~ 0.8). AGNs are known to have large variability am- 
plitudes in the X-ray band. For most objects in our sample the 
X-ray fluxes were determined based on one observation of the 
object which could well be away from its average flux state. 
Furthermore, different objects in the sample were observed 
with different instruments and telescopes having different ef- 
fective bandpasses, thus enlarging further the scatter between 
the objects. For most objects, the X-ray luminosity measure- 
ment is not contemporaneous with the time lag measurement 
(also during a monitoring campaign of months and years the 
X-ray luminosity could change over a large range) and this 
can add to the scatter in the relation. Nevertheless, the larger 
slope for the X-ray luminosity may be real, since it is in ac- 
cordance with the recent result that a ox in AGNs (a measure 
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certainties represent the standard deviations of the different 
results. We note that for the X-ray luminosity the two meth- 
ods give somewhat different slopes. The detailed results for 
each method are given in Tableland the preferred choice of 
slope and normalization should depend on the application it is 
used for. Simply averaging the measurements obtained from 
the BCES and FITEXY methods for the relation between the 
BLR size and the optical luminosity at 5 100 A, using one data 
point per object (sixth line of Table[5), we find: 



Rblr 
10 It days 



(2.23 ±0.21) 



( XL x {5\Wk) 
\ 10 ;i 



ergss 



0.69±0.05 



(2) 



When using all individual measurements for each object (fifth 
line of Table[3} the mean slope for the two methods is 0.65 ± 
0.04 and the mean normalization is 2.55 ±0.13. We also find 
that the fits indicate an intrinsic scatter of ~ 40% in these 
relations. This is the first time this intrinsic scatter has been 
quantified. 

Our results reflect on the naive theoretical predicted slope 
for the relation between the BLR size and the luminosity of 



a = 0.5 which is based on the assumption that all AGNs have 
the same ionization parameter, BLR density, column density, 
and ionizing spectral energy distribution (SED). The fact that 
for most energy bands the slope is different from a = 0.5 indi- 
cates that at least for some of these characteristics the simple 
assumption is not valid and they probably show an evolution 
along the luminosity scale. 
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